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Abstract 

Thermodynamic and structural properties of the tangent diatomic fluid are studied in the frame- 
work provided by the Reference Interaction Site Model (RISM) theory, coupled with a Modified 
Hypernetted Chain closure. The enforcement of the internal thermodynamic consistency of the 
theory is described in detail. The results we obtain almost quantitatively agree with available 
or newly generated simulation data. We envisage the possibility to extend the consistent RISM 
formalism to generic, more realistic molecular fluids. 
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Tangent diatomics, constituted by two identical hard spheres whose centre-to-centre distance is 
equal to the diameter cr of each sphere, provide a relatively simple prototype model for molecular 
systems. They can be considered as a member of two more general classes, one constituted by 
hard dumbbells (or fused hard spheres), where each sphere composing the molecule can have 
a different cr diameter and the bond length (elongation) is generally a fraction of cr, and the 
other one constituted by two or more freely jointed hard spheres, a model widely used as a basic 
representation for chain-like molecules. 

The structural and thermodynamic properties of the tangent diatomic fluid including its phase 
behavior, have been widely analyzed until recently in terms of both computer simulations and 
liquid state theories. Early Monte Carlo calculations were carried out by Freasier et al. @]- 
In Ref. [3] Tildesley and Streett have used the Monte Carlo pressure to evaluate the constants 
entering an empirical analytic expression that accurately fits the compressibility factor as a func- 
tion of the density. More recent simulations can be found in Refs. B4dL^IlZD. The fluid-solid 
equilibrium has been investigated via Monte Carlo simulations and free energy calculations in 
Refs. flU. 

Theoretical studies have involved the Reference Interaction Site Model (RISM) theory of 



molecular fluids developed by Chandler and Andersen yjj. In particular, in Refs. lid 1 1 111 the 



tangent diatomic fluid has been analyzed in the context of the Chandler-Silbey-Ladanyi "dia 



grammatically proper" formulation of RISM (CSL, 11211 ). Calculations based on the contracted 
formalism known as "polymer-RISM" are reported in Ref. [7]. Among other theoretical studies 
we mention a closed form analytic theory for the structural functions 1 13], an Ornstein-Zernike- 
type integral equation theory with a non-spherical bridge function lll4ll . a modification of the 
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Verlet theory for hard spheres [15], a scheme based on the Born-Green- Yvon integral equation 
theory iflill . Th e eq uation of state of tangent diatomics has been investigated via different for- 
malisms (see pi 1 171 1 1811 and references therein). 



The aim of this Letter is to present the first (to the best of our knowledge) thermodynamically 
consistent study of the tangent diatomic fluid in the framework provided by the RISM theory. 
We document the accuracy of thermodynamic and structural properties obtained through such a 
scheme and envisage further applications to more realistic molecular fluids. 

The RISM formalism ifioll is a matrix generalization of the Ornstein-Zernike equation of sim- 
ple fluids 01911 . relating the set of site-site pair distribution functions to the corresponding direct 
correlation functions. The molecular geometry enters the theory through a matrix of intramolec- 
ular correlations that takes into account the rigid bonds among the various interaction sites of the 
molecule (see Ref. 12011 for a review of the method and applications). Because of the symmetry 
of the tangent diatomic molecule all site-site correlations are equal, so that the RISM equation 
assumes the simple form in k space: 

h(k) = [w(k) + l] 2 c(k) + 2p[w(k) + l]c(k)h(k) . (1) 

In Eq. (0} p is the molecular number density, h(r) = g(r) - 1 and c(r) are respectively the pair 
and direct correlation functions among any of the two sites of different molecules, and g(r) is the 
site-site radial distribution function. The function \v(k), that takes into account the intramolecular 
correlations, is written as: w(k) = sin{kL)/kL, where L = cr is the bond distance among the two 
spheres composing the model. 

The RISM equation must be complemented by a closure relation, that is usually taken to be 
the Percus-Yevick (PY) or the Hypernetted Chain (HNC) Jiol Eltl . Going beyond these basic 



approximations, we suggest here to adopt a Modified Hypernetted Chain (MHNC) closure 12211 : 
specifically, we assume that the well known exact expression for the radial distribution function 
of an atomic fluid lfl9ll . 

g(r) = exp[-/Mr) + W) ~ c{r) + E(r)] , (2) 

may be employed for the site-site structural functions of a molecular fluid. In the equation above 
v(r) is the interparticle potential, p 1 - l/k^T (where k% is the Boltzmann constant and T the 
temperature) and E{r) is the "bridge function". The quotation marks are used in this context 
because the theoretical framework for a rigorous definition of the bridge diagrams is provided by 



the CSL proper formulation of RISM 111 20 (see also the recent developments in Refs. 11 ll 12310 . 



As in the original MHNC scheme for atomic fluids 12211 . we then approximate E(r) by Ens(f), 
namely the bridge function of a hard sphere fluid of some effective packing fraction /7ns , in the 
parametrization of simulation data provided by Verlet and Weis [24]. 

In the MHNC scheme for atomic fluids, 77hs is adjusted to enforce the thermodynamic con- 
sistency of the theory, as for instance by requiring the equality between the virial and the com- 
pressibility equations of state. For a generic molecular fluid the virial equation of state cannot be 
deduced in terms of site-site radial distribution functions 12511 : however, for hard dumbbells, and 
hence for the model at issue, an expression for the excess free energy per particle has been de- 
rived by Lowden and Chandler j26ll in terms of an integral over cr of the value of g(r) at contact, 
S(cr + ): 

/SA " = 8/rp f" d<r'(<r'fg(r = cr' + ;p, cr') . (3) 
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Figure 1 : MHNC, HNC and PY compressibility factor of the tangent diatomic fluid. Monte Carlo data are taken from 
Ref. d]. 



This expression for A ex leads, upon derivation with respect to the density, to a "virial-like" equa- 
tion of state, that can be made equal (with a suitable choice of ^hs) to the pressure calculated 
from the compressibility route: 

1 dp 

--£ = l +ph (k = Q). (4) 

The coupled equations (Q]) and (O have been solved through standard numerical methods. We 
employ a discrete grid of 8192 points, with a r spacing Ar = 0.005cr. Complementary calcu- 
lations with 2 15 points and Ar = 0.002cr show no appreciable difference with the smaller grid. 
In order to enforce the thermodynamic consistency, the calculations are repeated at each density 
for different values of 77ns till the pressures coming from Eq. © and (|4]i coincide within a ~ 2% 
numerical accuracy. 

We report in Fig. [T|the RISM/MHNC compressibility factor /3P/p for the tangent diatomic 
fluid, gauged against Monte Carlo data Q3J] , and in comparison with the HNC approximation, 
(corresponding to the assumption E{r) — in Eq. (fJJ) and the PY closure (in which one assumes 
c(r) = outside the hard-sphere core). As visible, the thermodynamically consistent MHNC 
reproduces quite well the Monte Carlo equation of state over practically the whole fluid den- 
sity range (the freezing threshold of the model being estimated at per 3 ~ 0.53 Jll]). The HNC 
virial-like and compressibility equations of state bracket the simulation data; as an example, they 
exhibit an inconsistency of about 20% of their average value at per 3 = 0.4. The HNC pressure 
obtained from a closed expression derived in Ref. [ 27] (not shown in the Figure) is close to, but 
systematically above, the virial-like equation results. Both PY routes lay systematically close 
one to each other and slightly underestimate the simulation results. In order to complete our 
picture, the comparison with other microscopic theories shows that the MHNC approach and the 
compressibility route in the Born-Green- Yvon (BGY) theory (see Figure 5 of Ref. [16]) share 
the same level of accuracy for the pressure up to per 3 ~ 0.4, whereas the BGY virial route un- 
derestimates systematically the simulation data from per 3 ~ 0.3 onwards. The CSL approach 
reproduces quite well the Monte Carlo equation of state (see Figure 3 of Ref. |@]), although it is 
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slightly less predictive than the MHNC in the high density fluid regime. 

The quality of our results can also be appraised from Figure [2] where we display the MHNC 
input functions for the consistency procedure, namely the excess free energy calculated from 
Eq. ([3]) and the k — > limit of the static structure factor, S (0) which constitutes the right hand 
side of Eq. (@). The MHNC results are presented along with the HNC and PY predictions, and 
compared with the corresponding analytical functions derived in Ref. H from an accurate fit of 
simulation data for the pressure. We observe that, although all closures give a relatively accurate 
reproduction of simulation data, the integration of the compressibility and the derivation of the 
free energy eventually give rise to the relative spread of theoretical predictions for the pressure 
already shown in Fig.Q] 

MHNC and PY predictions for the site-site radial distribution function g(r) are displayed in 
Figure [3] and compared with Monte Carlo data generated in this work. Standard Monte Carlo 
simulations are generally carried out on samples composed of N - 500 molecules enclosed in a 
cubic box with periodic boundary conditions. We have checked in these conditions the negligible 
influence of the box size and of the spatial mesh of structural functions through several runs with 
4000 particles and by using spatial grids as fine as O.OOlcr spacing. It appears from Figure[3]that 
in the intermediate-to-high density regime (per 3 > 0.3) the MHNC correlations are practically 
superimposed to the simulation results. Only at the highest density investigated, per 3 = 0.5, the 
"exact" g(r) displays a shoulder before the cusp at r — 2cr preluding to the freezing of the fluid, 
whereas the MHNC is less sensitive to this emerging feature. At per 3 = 0.1/0.2, the theoretical 
correlations agree with simulation data for distances around the cusp onwards, whereas they 
appear slightly distorted in the region cr < r < 2cr. This artifact is produced by too high a value 
of 77HS, necessary to impose the thermodynamic consistency, As a result, the contact value of 
g(r) is slightly underestimated in the dilute fluid regime. The overall quite good reproduction 
of g(r) at contact shown in Fig. [3] is noteworthy (the discrepancy with simulation data barely 
exceeding ~ 5% in the worst case investigated, i.e. at per 3 = 0.1), since g(cr + ) constitutes an 
important parameter entering not only Eq. (|3), but also other approximate expressions for the 
equation of state of tangent diatomics [16]. As visible from the figure, the MHNC and the 
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Figure 3: MHNC (full lines), PY (dashed lines) and Monte Carlo (symbols) site-site radial distribution functions for the 
tangent diatomic fluid at various densities. Theoretical predictions lay close to the corresponding simulation results and 
are superimposed at per 3 =0.3. 



PY predictions are indistinguishable at per 3 = 0.3. The PY is sligthly less accurate than the 
MHNC outside this density range, and in particular, as already observed in Ref. fll, tends to 
overestimate (underestimate) the contact value of g{r) at lower (higher) densities. In comparison 
with our predictions, results from BGY Biol and CSL |6] theories show a better agreement with 
simulations data only in the low density regime. 

The body of theoretical results presented in this Letter appears of remarkable accuracy, espe- 
cially in consideration of the relative simplicity of the scheme adopted. The thermodynamically 
consistent MHNC improves on previous RISM results and exhibits an accuracy similar to the 
more rigorous but comparatively more complex CSL formalism. We have also several evidences 
testifying that the MHNC approach positively predicts the thermodynamic and structural proper- 
ties of hard dumbbells for generic elongations other than L = <r B28I1 . The present study paves the 
way for further applications to more realistic systems. In fact, although a consistency procedure 
for a generic molecular model cannot hinge on Eq. ©, since this relationship holds strictly for 
the hard-dumbbell fluid, other consistency schemes, based e.g. on the "energy route" to ther- 
modynamics, can be implemented in such cases. Investigations on these topics are currently 
underway. 
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